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(54) IMethod and apparatus for fast aerial image simulation 



(57) The present invention provides a method and 
apparatus for simulating an aerial image projected from 
an optical system, wherein the optical system includes 
a pupil and a mask. In general, the method comprises 
the steps of obtaining parameters for the optical system, 
calculating a kernel based on an orthogonal pupil pro- 
jection of the parameters of the optical system onto a 
basis set, obtaining parameters of the mask, calculating 
a vector based on an orthogonal mask projection of the 
parameters of the mask onto a basis set, cateulating a 
field Intensity distribution using the kernel and the vec- 
tor, and obtaining aerial image data from the field inten- 
sity distribution. 
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Description 

BACKGROUND OF THE INVENTION 

5 [0001] In lithography, an exposure energy, such as ultraviolet light that is generated in an optical system, is passed 
from an aperture of the systenn through a mask (or reticle) and onto a target such as a silicon substrate. The mask 
typically may contain opaque and transparent regions fomned in a predetermined pattem. The exposure energy exposes 
the masl< pattem, thereby fomiing an aerial image of the mask. The aerial image is then used to fomn an image onto 
a layer of resist formed on the target. The resist is then developed for removing either the exposed portions of resist 

10 for a positive resist or the unexposed portions of resist for a negative resist. This forms a patterned substrate. A mask 
typically may comprise a transparent plate such as fused silica having opaque (chrome) elements on the plate used 
to define a pattern. A radiation source illuminates the mask according to well-known methods. The radiation transmitted 
through the mask and exposure tool projection optbs fonns a diffraction-iimited latent image of the mask features on 
the photoresist. The patterned substrate can then be used in subsequent fabrication processes. In semiconductor 

IS manufacturing, such a patterned substrate can be used In deposition, etching, or Ion implantation processes, to fonm 
integrated circuits having very small features. 

[0002] In a manufacturing process using a lithographic projection apparatus, a pattem In a riiask is Imaged onto a 
substrate, which is at least partially covered by a layer of radiation-sensitive material (resist). Generally, lithographic 
patterning processes are understood by those who practice the profession. Information regarding exemplary processes 
^ may be obtained, for example, from the book "Microch^ Fabrtoation: A Practical Guide to Semiconductor ProoesslngT, 
Third Edition, by Peter van Zant, McGraw Hill Publishing Co.. 1 997 ISBN 0-07-067260-4. 

[0003] As the size of lithographically fabricated structures decreases, and the density of the structures increases, 
the cost and complexity of designing masks additionally increases. One method of reducing costs of lithographic fab- 
rication is by optimizing the lithographs design with a lithographic simulation step prior to the actual manufacturing 

ss step. One specific method of lithographic simulation Is drawn to simulating the aerial image of the mask. The aerial 
Image Is defined as an intensity distribution of light just prior to reaching the resist on a surface of a substrate, when 
the substrate is exposed via the mask i.i an exposure apparatus. In order to simulate the aerial image, a layout of a 
mask and exposure conditions (non-limiting examples include NA: Numerical Aperture, a (sigma): Partial Coherence 
Factor) of the lithographk^ apparatus are typically required as input parameters. 

50 [0004] Lithographic apparatus may employ various types of projection radiation, non-limiting exEimples of whteh in- 
clude light, ultra-violet ("UV") radiation (including extreme UV ("EUV"), deep UV ("DUV"). and vacuum UV ("VUNT)). 
X-rays, ion beams or electron beanrs. The following have been considered exemplary exposure sources. Light may 
generally refer to certain mercury emissions, i.e., wavelengths of 660 nm for the f-ilne, 436 nm for the g-llne, and 405 
nm for the h-iine. Near-UV or UV generally typteally refer to other mercury emissions, i.e., 365 nm for the i-line. DUV 

55 generally refers to excimer laser emissions, such as KrF (248nm) and ArF (1 93nm). VUV may refer to excinner laser 
Fg, i.e.. 167nm, Arg, I.e., 126nm, etc. EUV may refer to 10-15 nm. This last portion of the electromagnetic spectmm is 
very close to "soft X-rays" but has been named as "EUV", possibly to avoid the bad reputation of X-ray patteming. Soft 
X-rays may refer to 1-15 nm, whteh may typically be used in X-ray lithography. 

[0005] Depending on the type of radiation used and the particular design requirements of the apparatus, the projection 
40 system may be for example, refractive, reflective or catadioptrto, and may comprise vitreous components, grazing- 
Incldence min-ors, selective multi-layer coatings, magnetic and/or electrostatic field lenses, etc; for simpltelty, such 
components may be loosely referred to In this text, either singly or collectively, as a "lens". 

[0006] When the resist is exposed by the aerial image, there is an additional variable In that some of the exposure 
light Is reflected back by the surface of the substrate, and then absort3ed by the resist. Accordingly, not only the resist 

45 characteristtes (regarding, for example, refractive Index: Dlirs A,B,C) but also parameters regarding the characteristics 
of the substrate (e.g., refractive Index) should be Included as input parameters for simulating the latent Image. 
[0007] The so-called Hopkins model treats the electric field fomiing the image typically as a scalar and assumes the 
object being imaged is sufficiently thin so that its effect on the incident field is represented by a multiplicative function. 
It is advantageous to perfomr) the image fomnatfon analysis in the Fourier domain (frequency space) In order to deal 

so with the pupil function of the Imaging system rather than the amplitude response function, and with the angular distri- 
bution or "effective source" rather than with the mutual intensity. 

[0008] There are several computer programs commercially available that calculate aerial images based on the Hop- 
kins model. For example, the University of California at Bericeiey. Department of Electrical Engineering and Computer 
Science, Bericeley, Calif., 94720, offers a program called SPLAT. 
55 [0009] The Hopkins model is usedto model the imaging of drawn design features under partially coherent Illumination. 
A major problem in modeling aerial images under partially coherent Illumination is the necessity to superimpose and 
add the effect of each Individual illumination source that makes up the partially coherent source. In the Hopkins model, 
a two-dimensional by two-dimensional transmission cross-coefficient function ("TCC") is pre-calculated, which captures 
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all the effects of the lithographic projection apparatus, incfuding NA. sigma. etc. As taught, for example in Bom & Wolf 
p. 603. once a TCC Is l<nown. systems with partially coherent Illumination can be modeled by integrating the TCC over 
the Founer transfomi of the transmission function for the geometrical layout feature under illumination 
[0010] Fundamentally, the TCC Is a two^limensional by two-dimensional correlation function with a continuous set 
s of arguments. In practice, an assumption can be made that the feature patterns to be imaged are periodic in space 
For such penodic patterns, the TCC has a large, but discrete, set of arguments. The TCC can then be represented as 
a matrix, with discrete columns and rows. For typical features interesting to the lithographer the size of this matrix Is 
tremendous and restricts tlie scale and size of features that can be simulated. It Is a purpose of any simulation algorithm 
10 of accural ^ '° '^'""^ °^ ^"'^ approximation, while retaining a reasonable degree 

An exemplary projection lithography system is illustrated in FIGS. 1 A and 1 B. In FIQ. 1 A, light from illumination 
source 102 is focussed by condenser 104. The condensed light passes through the mask 106. then through the pupil 
108 and oiito the substrate 110. As shown in FIG. IB, substrate 110 may comprise a top anti-reflective coating 112. 
a resist 118 a bottom anti-reflective coating 117. a top substrate layer 114 and a plurality of other substrate layers. 
As Indicated In FIQ. 1 B, the focal plane may tie within the resist 118. 

I0012J In the past, several numerical techniques have been applied to reduce the size of the TCC to reasonable 
scal^. in one Instance singular value decomposition has been applied to decompose the TCC Into its eigenspectrum 
sort the resulting eigenvectors In decreasing magnitude of their eigenvalues, and only retain a finite number ofeigen- 
vectora in order to approximate the TCC. An exemplary method for optical simulation for the system of FIQ. 1 A that 
uses the Hopkins model, Is illustrated in FIQ. 2. As illustrated In RQ. 2. at step S202. the lithographic projection ap- 
paratus and mask parameters are input fn the system. At step S204. the TCC is approximated by an eigenvector 
diagonalizatlon. At step S206, a two-dimensional table of convolution integrals of the eigenfunctions with a discrete 
set of normalized rectangles is oateulated and tabulated. The article, "Fast, Low-Complexity Mask Design", by N Cobb 
K t J 3^3-326. the entire disclosure of whfch is incorporated by reference, teaches an exemplar/ 

method for completing step S206. Since the resulting table is on a discrete grid of possible layout-features, to attain 
an acceptable resolution, the table must be large. At step S208. since the TCC has been approximated, the Illumination 
system can be modeled, i.e. the aerial image can be simulated by combining each two-dimensional convolution Integral 
conesponding to each respective rectangle in a given proximity window. 

^013J This method has two drawbacks. First, an eigenvector diagonalizatlon is a numerically expensive operation 
TTius, for a user to change lithographte projection apparatus conditions, such as NA. sigma, illumination type or lens- 
aberration an expensive recateulation of the TCC approximation that limits the use of the simulation tool is required 
Second, the eigenspectnim of the TCC is a two-dimensional function. Representing the associated field vectors of the 
geometncal pattem under illumination requires a two-dimensional lookup table that limits the speed otthe Held calcu- 
ation m the Hopkins algorithm. Specifically, since the field vectors are a two-dimensional table, two pointers are needed 
to access any respective address of data, whteh increases access time. Further, since the field vectors are a two- 
dimensional table, they are stored in a portion of the cache. Using a look-up table from a DRAiW in a present day CPU 
IS more experwive than perfomilng multiplication because the cache miss penalty is larger than the time to perform 
the multiplication. 

100141 Accordingly there is needed a method and apparatus for simulating a projection lithography system using 
the Hopkins model, that does not use a numerically expensive eigenvector diagonalizatlon operation. What is further 
needed is a method and apparatus for simulating a projection lithography system using the Hopkins model that permits 
a user to change lithographic projection apparatus condittons. such as NA. sigma, illumination type or lens-aberration 
without requiring an expensive recalculation of the TCC approximation. 

BRIEF DESCRIPTION OF THE INVENTION 

[001S] II is an object of the present invention to provide a method and apparatus for simulating a projection lithography 
fno^T "f """"""^ ^^^^ "^^ ^ numerfcally expensive eigenvector diagonalizatlon operation 

30 mh«^ K another object of the present invention to provide a method and apparatus for simulating a projection 
litl^ography system usingthe Hopkins model that pemiits a usertochange lithographic projection apparatus conditions 
without requlnng an expensive recalculation of the TCC approximation. 

I^IV .'I'®.**'" object of the present invention to provide a method and apparatus for predicting the intensity 

field distnbutlon faerial image") at the surface and thmughout the resist film of a substrate being Irradiated ir^^aZ 
graphic imaging process. 

55 [0018] In accordance with the foregoing objectives, the present Invention writes the TCC as T(q'q) where q and a' 
LTximlZ^^^ frequencies In two-dimensional Founer space, q=(q,. q,), q-(q. . q-^). The present Invention ap- 
proximates he TCC e)rtremely well as a bilinear form of a basis function with kernel A[y]. Furthemiore, the present 
Invention utilizes a set of orthogonal polynomials. The kemel A[IJ] represents a small matrix that efficiently approximates 
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the TCC for a broad range of Illumination conditions used In modenn lithographic processing. The present inventbn 
uses a lower number of arithmetic operations as a result of the use of a one-dimensional look-up table as opposed to 
the two-dimensional look-up table in the prior art. Furthennore. the one-dimensional look-up table achieves a higher 
cache hit ratio than a two-dimensional look-up table of the same resolution, thereby providing a more efftelent system. 

5 [0019] In general, In one embodiment, the inventfon features a method of simulating an aerial image projected from 
an optical system, the optical system including a pupil and a mask plane, the method comprising the steps of providing 
a mask to the mask plane, obtaining parameters for the optical system, calculating a kemel based on an orthogonal 
pupil projection of the parameters of the optical system onto a basis set, obtaining parameters of the mask, calculating 
a vector based on an orthogonal mask projection of the parameters of the mask onto a basis set, calculating a field 

10 Intensity distribution using the kemel and the vector, and obtaining an aerial image from the field Intensity distribution. 
[0020] In one embodiment of the invention, the parameters for the optteal system include abenrations. 
[0021] In another embodiment of the invention, the step of calculating a kemel con-esponding to the parameters of 
the optical system includes the step of generating a table of indefinite integrals based on a recurrence over a seed 
array of incomplete gamma functions. 

IS [0022] In yet another embodiment of the invention, the step of calculating a kernel corresponding to the parameters 
of the optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points In the pupil of the optteal system, wherein the optical system Is In-focus. More particularly, 
It further comprises the step of combining sample weights of an Illuminator profile of the optical system with the an^ay 
of orthogonal pupil projection coefftoients. 

^ [0023] In still another embodiment of the invention, the step of calculating a kemel corresponding to the parameters 
of the optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in the pupil of the optical system, wherein the optfcal system Is either not in-focus or has 
aberrations. More particularly, it further comprises the step of combining sample weights of an illuminator profile of the 
optical system with the array of orthogonal pupil projection coeffk^lents. 

25 [0024] In still yet another embodiment of the invention, the step of calculating a kemel con-esponding to the param- 
eters of the optical system further includes the step of: tabulating an array of orthogonal pupil projection coefficients 
corresponding to respective points in the pupil of the optical system, wherein the optical system accounts for effects 
of diffusion of photoactive compounds in the resist. More particularly, it further comprises the step of combining sample 
weights of an illuminator profile of the opttoal system with the an^y of orthogonal pupil projection coefficients. 

30 [0025] In a further embodiment of the invention, the step of calculating a vector connesponding to the parameters of 
the mask further includes the step of specifying a proximity window within the mask for geometric sampling. 
[0026] In still a further embodiment of the invention, the step of calculating a vector corresponding to the parameters 
of the mask further includes the step of decomposing a geometric patlem of the mask into a disjoint set of rectangles 
and tabulating an array of projections of the rectangles within a proximity window. More particularty, the step of calcu- 

35 lating a vector conresponding to the parameters of the mask further Includes the step of conBCting the an^ay of projec- 
tions of the rectangles based on the type of mask. 

[0027] In general, in another aspect, the invention features a method for simulating an aerial image projected from 
an optical system, the optical system including a pupil and a mask plane, the method comprising the steps of: providing 
a mask to the mask plane, obtaining parameters for the optteal system, obtaining parameters of the mask, projecting 

•to the components of electric field vectors to an orthogonal polynomial basis as described in further detail below, the 
polynomials including polynomials associated with the parameters of the optical system, and approximating the trans- 
mission cross-con-eiation function associated with the optical system based on the orthogonal projection of polynomials. 
[0028] In general, in still another aspect, the inventton features a simulation device operable to simulate an aerial 
Image projected from an optteal system, the optteal system Including a pupil and a mask plane, the simulation device 

45 comprising a first parameter obtalner for obtaining parameters of the optical system, a first calculator for calculating a 
kemel based on an orthogonal pupil projection of the parameters of the optical system onto a basis set, a second 
parameter oblainer for obtaining parameters of a mask provided to the mask plane, a second calculator for calculating 
a vector based on an orthogonal mask projection of the parameters of the mask onto a basis set, a third calculator for 
calculating a field intensity distribution using the kemel and the vector, and an aerial Image obtalner for obtaining aerial 

50 image data from the field intensity distribution. 

[0029] In one embodiment of the Invention, the first, second, and third calculators are the same calculator. 
[0030] In another embodiment of the invention, the parameters for the optical system include aberrations. 
[0031] In yet another embodiment of the Invention, the first calculator is operable to generate a table of indefinite 
Integrals based on a fast recurrence over a seed an^ay of incomplete gamma functions. 

55 [0032] In still another embodiment of the invention, the first calculator is operable to tabulate an array of orthogonal 
pupil projection coefficients corresponding to respective points in the pupil of the optical system, wherein the optical 
system Is in-focus. 

[0033] Instill yet another embodiment of the invention, thefirst calculator Is operable to tabulate an array of orthogonal 
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pupil projection coefficients corresponding to respective points In the pupil of the optical system, wherein the optical 
system is either not in-focus or has aberrations. 

[0034] In a further embodiment of the Invention, the first calculator is operable to tabulate an array of orthogonal 
pupil projection coefficients corresponding to respective points in the pupil of the optical system, wherein the optical 
system accounts for effects of photoactive compound diffusion In the resist. 

[0035] In still a further embodiment of the Invention, the second calculator is operable to specify a proximity window 
wUhin the mask for geometric sampling. 

[0036] In yet a further embodiment of the invention, the second calculator Is operable to decompose a geometric 
pattem of the mask into a disjoint set of rectangles and tabulate an an^y of projections of the rectangles within a 
proximity window. More particulariy the second calculator is operable to con-ect the array of projections of the rectangles 
based on the type of mask. 

[0037] In general, in yet another aspect, the invention features a simulation device operable to simulate an aerial 
image projected from an optical system, the optical system including a pupil and a mask plane, the simulation device 
comprising a first parameter obtainer for obtaining parameters of the optical system, a second parameter obtalner for 
obtaining parameters of a mask provided to the mask plane, a first calculator for orthogonally projecting polynomials, 
the polynomials Including polynomials associated with the parameters of the optical system, and a second calculator 
for approximating the transmission cross-coefficient associated with the optical system based on the orthogonal pro- 
jection of polynomials. 

[0038] Additional advantages of the present Invention will become apparent to those skilled in the art from the fol- 
lowing detailed description of exemplary embodiments of the present Invention. The invention Itself, together with 
further objects and advantages, can be better understood by reference to the following detailed description and the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0039] The accompanying drawings, which are incorporated In and form a part of the specification, illustrate embod- 
iments of the present invention and., together with the description, serve to explain the principles of the invention. 

In the drawings: 

[0040] FIGS. 1 A and 1 B illustrate a conventional optical lithographic system. 

[0041] FIG. 2 Illustrates a conventional method for simulating an aerial image projected from an optical system. 
[0042] FIG. 3 illustrates a conventional method for dividing a pupil area Into a discrete areas. 
[0043] FIG. 4 illustrates the method for simulating an aerial image projected from an optical system in accordance 
with an exemplary embodiment of the present invention. 

[0044] FIG. 5 Illustrates an exemplary method for generating a kernel in accordance with the present invention. 

[O045] FIG. 6 Illustrates an exemplary method for generating a vector in accordance with the present invention. 

[0046] FIG. 7 illustrates an exemplary method for generating a table of indefinite integrals based on a fast recun-ence 

over a seed array of incomplete gamma functions In accordance with the present invention. 

[0047] FIG. 8 illustrates an exemplary mask plane used in conjunction with the present invention. 

[0048] FIG. 9 Illustrates an exemplary method for calculating an Intensity profile and intensity gradient In accordance 

with the present Invention. 

[0049] FIG. 1 0 illustrates an efficient algorithm in accordance with the present invention. 

[0050] FIG. 1 1 A illustrates the coefficient matrix A^ in accordance with the present invention. FIG. 11 B illustrates an 
exemplary representation of the coefficient matrix A^ In accordance with the present Invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0051] In the following description, for the purposes of explanation, numerous specific details are set forth in ortler 
to provide a more thorough understanding of the present Invention. It will be apparent, however, to one skilled In the 
art that the present invention may be practiced without these specific details. 

[0052] The present invention circumvents the above-described problems of the prior art systems by applying a fun- 
damentally different technique for the approximation of the TCC. Specifically, the present Invention closely and effi- 
ciently approximates the TCC as a bilinear fomi of a basis function with kemel A[lj]. 

[0053] A general oven/lew of the operation of the simulation system In accordance with the exemplary embodiment 
of the present invention is illustrated in FIG. 4. As Illustrated in FIG. 4 of the present invention, in step S402, the 
lithographic projection apparatus parameters are entered. In step S404, a kernel A[lj] is calculated. In step S406, the 
mask parameters are entered. In step S408, vectors corresponding to mask parameters are calculated. In step S410 
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the intensity profile of the pupil Is calculated using the kernel A[ijl and the vectors corresponding to the mask parameters. 
In step S412, a simulated aerial image of the mask is produced. 

[0054] Prior to discussing the operation of a simulation system in accordance with the exemplary embodiment of the 
present Invention in further detail, a brief but concise discussion of the mathematical underpinnings will be provided. 
5 Specifkjally, the spectral basis, including its analytical representation, its basic properties, and Its appfeatlon with re- 
gards to an indefinite integral will be discussed. Further, the Fourier transfonn ("FT*) of the basis, including Its analytical 
representation, and its application with regards to an indefinite integral will be discussed. Still further, miscellaneous 
provisions including orthogonality under a FT will be discussed. 

10 SPECTRAL BASIS 

Analytteal Representation 

[0055] A mathematical function provides a mapping of a continuous set of input values to a continuous set of output 
IS values. This mapping may be generated In a number of ways. One way is the use of an analytical expression, such 
as sin(x), or log(x). Another way is by tabulation and interpolation of intemiediate values. Yet another way is to represent 
a function in temris of a linear superposition of a set of known basis functtons. For an orthogonal set of basis func^ns 
such a representation is called a spectral representation. 

[0056] The spectral basis considered in the following discussion is represented as: 

20 

<i1[a]"> = (a/^/^i){-i)^'' exp[7^ / 2]d%%xf['f] (1) 

25 which may be factored into one-dimensional contributions. In equation (1), r = a (r-rQ) for sampling at a point Tq, and 
Xg = a,( -i- ay The one-dimensional basis along the x-dimension follows as: 
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Basic Properties 
[0057] 

1 . Recurrence 

2. Derivative 

^"^(al"^" § aXnfa.^/|l^-<n[a^1.jl'> W 

3. Normalization 

The basis is orthogonal with normalization constants: 

(7. (r)- 7. (r)) = with a=(ax, a^) (5) 



55 Indefinite integral 

[0058] The Indefinite integral overthe spectral basis may be related to an incomplete Gamma function. For Manhat- 
tan-style layout designs, the integration boundaries are along the Cartesian coordinate axes and pennit a factorization 
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into independent one-dimensional integrations. The following definition of the incomplete Gamma-function v(ax) 
used: *^ ' ' 
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[0059] Recurrence 



Let x= ax and 



(6) 



(7) 
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[00601 Note that J(a;,,„](-i) = (-1 * ' J^^^ (i). The above Integral may then be evaluated 

•[<i,(>fo.^l)=J[a.oi<^l)-'taO)(*^0) 



(8) 



as 



(10) 
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[0061 ] The Fourier transfomrj of the spectral basis may be expressed as: 
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Indefinite Integ ral 



(11) 



(12) 



[0062] Similar to equation (3). the indefinite Integral over the FT of the spectral basis may also be related to ... 
moomplete Gamma function. Due to the circular aperture of the pupil however the two-dimensional integration cannot 
be simply factorized Into one-dimensional products. The present invention uses a hybrid approach, whereby the one- 
d^ensional integration along one axis, for example the x-axis, is perfonned analytically similar to equation (3) On the 
other hand, the integration along the other axis, for example the y-axis. is performed numerically. The following 
siders the one-dimensional analytical integration of the FT of the spectral basis. Let 
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:= (13) 

a 

= a J*'d?(7^.^k)-a |*''rfg{7^,^w) (15) 



Recurrence 

15 

[0063] Let /c= /c/aand 



^,*,.«3W := [-^]^[0J„' ArexptV2]arexp[-f^] (16) 



then the above integral may be evaluated as 

25 



A recurrence for J\a^m] (k) follows via 
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- J'^^(0"exp[r'/2]a,arexp[-/^l (18) 



35 

as: 



40 

[0064] This relation is identical to equation (1 0), with p = 1 /a and an additional factor of (-/) per recurrence step, 
initialization 

45 [0065] The above seed-elements c/fjo are similar to the previous section and may be evaluated similar to equation 
(11): 



t2 

jf 1 2 ' ar .m+i 

so 



[0066] With the above refreshing review, the operation of a simulation system in accordance with the exemplary 
embodiment of the present invention will now be provided in further detail. Table 1 , below, provides notations for symbols 
S5 and their respective descriptions that are used throughout the description of the present invention. 
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Symbol 



(Ky) 



m,l 



M 
S 



Description 



2-D position vector fixed z-plane 
2- O wave vector fixed z-plane 
plane wave entering mask 

scalar field at mask exit plane at nonnal Incidence, d = 0 
scalar field at pupil entrance plane, general Incidence angle 
scalar field at pupil exit plane, general Incidence angle 
scalar field at wafer entrance plane, general Incidence angle 
Mask transfer matrix 
Pupil transfer matrix 
Magnifksation factor (0.2, 0.25, etc) 
Scaled coordinate system at wafer plane 
Scaled coordinate system at pupil plane 
Partial Coherence 



[0067] The following discussion uses DIrac notation for a basis^independent representation of field vectors. In this 
20 notation, scalar fields in their dual spaces are written as <<|>l. I(|»>. The following are useful identities: 



25 



30 



(21) 
(22) 
(23) 
(24) 



35 



(hk) = (h'k) 



(25) 



(rik-k') = <k'lr)<rlk> (26) 
[0068] Completeness of an orthogonal nomnaiized basis ^ Implies 

/ d^l^X^ 1 = 1 (27) 

[0069] Spectral expansions used in the following satisfy the above completeness relation. 

Lens Model 

[0070] When modeling a lens, monochromatic illumination at wavelength X is assumed. In the following the refractive 
Index in vacuum 7)/= 1 is assumed. The wave-vector may be defined as: 



40 



45 
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I, 27C 

[0071] Further, the apertures may be defined as (w= substrate, s = Source): 



(28) 
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NA^ = Pf sin ^ sma^ (29) 



NAg - Hf sin - sin (30) 
NA^ 



10 



[0072] The maximunn wave-vector transmitted by a pupil of aperture NA thus follows as 

k^=d^sma^=d^ NA (32) 
[0073] A magnification factor M may be defined as the ratio of the bacl<- vs. the front-focal length of the objective lens: 

Af:=J (33) 

[0074] Where dq and d^ cooespond to the focal lengths identified In Fig. 1 A. The substrate plane may be adjusted 
such that It lies in the focal plane of the lens, thereby providing: 

/ do 

[0075] The mathematical model is expressed in temns of coordinates kg, r, k, The numerical model however 
utilizes re-scaled coordinates indicated by "har notation. The scaled coordinates relate to the normalized coordinates 
via 



15 



20 



25 



35 



40 



45 



K '^^'^^m Mask plane (35) 



:= l<o sin d = Pupil plane (36) 



Jfc^'' koAM (37) 



so 



A kfl sin 9- 

s •= - -4 = - Source plane (38) 

The current lens model in accordance with the present invention assumes = do. 
DIffractive Elements and Lenses 

[0076] The propagation of light from the exit plane of a diffraction element, such as for example a mask. Ir) to a 
subsequent input plane <<t>' I /c> at distance d^ may be expressed in terms of the free-space Greens function of the 
Helmholtz equation < £ I G I r*)(r = {x,y.z) in three dimensions) as 
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(39) 



with 



15 



20 



</1Glf'> = 



4n\r't\ 



(40) 



10 [0077] Let z and t be on fixed planes separated by a distance d/ and r = (x.y) a two-dimensional position vector on 
each respective plane. A two dimensional approximation follows via 



k^lr-/;'l = k^7df+(r-0^ (41) 
[0078] Assume a fer-fleld approximation Ir - r' I d, and approximate by keeping to first order of a Taylor expansion 



1 + 



(42) 



25 



30 



[0079] The first term yields an irrelevant constant phase shift and will be ignored in the following. The two-dimensional 
field propagator in the far-field regime then follows as 



<rlG ir') = 



4ndt 



(43) 



Lens Element 



[0080] The propagation of light from the input plane to the exit plane of a lens element may be expressed In terms 
of the lens propagator L as 



40 



<rl<|»''> = /tyr'<rlLlr*><r'l<J''*>=<rlLI<^'''> 
where the lens-propagator In the paraxial approximation may be expressed as: 



(44) 



45 
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(46) 



[0081] Returning back to FIG. 4, the simulation method of the present invention may be implemented using any 
processor or dedicated circuit that is operable to manipulate data. Non-limiting examples include programmable logic 
an^ays, microprocessors used in conjunction with a memory, or computers. Step S402 may be perfonned by any device 
operable to Input lithographic projection apparatus parameters. Non-limiting examples Include stepper/scanners, mask 
Inspection tools, and external memories containing the appropriate data. Step S406 may be performed by any device 
operable to input the mask parameters. Similar to step 8402, non-limiting examples Include stepper/scanners, mask 
Inspection tools, and external memories containing the appropriate data. Further, a single device may be used to input 
both the lithographic projection apparatus parameters of step S402 and the mask parameters of step S406. Still further, 
steps S402 and S406 may be perfonned at the same time. Similarly, steps S404 and 8408 may be performed at the 
same time. 

[0082] Step 8404 will now be further described with respect to FIG. 5. At step 8502, an area of the pupil of the 
illuminalor is sampled and broken up into a discrete grid. At step 8504, the illumination profile of each discrete point 
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corresponding to the sample area is determined, and given a numerical weight according to the relative intensity of 
the source at the sampled location. At step S506, an anBy of seed-elements for the semi-Indefinite integrals calculated 
In equation (1 6) of each discrete point within a pupil grid con-esponding to the sample area Is tabulated. Each seed- 
element is an incomplete gamma function J^o.m]. such as is described above with respect to equation (1 9) and equation 

s (20). At step S508, the pupil model is detemiined for the simulation method. If the pupil model is for a zero defocus 
pupiL i.e.. for a lens having no aberration, Is perfectly in focus, and/or has no phase distortion in the pupil plane, then 
step SSIOzero be performed. If the pupil model has optical path differences, has aberrations, is not perfectly in 
focus, and/or has phase distortion in the pupil plane, then step SSIOopd "^^V ^® perfomied. Further, if the pupil model 
accounts for diffusion of the photoactive compound in the resist after Image exposure, then step SSIOdttr may be 

10 perfonned. In each one of steps SSIO^ero' ^^^^opd' ^^^^tm* ^^^^ shifted by S/, of pjisjt Is calculated, 
wherein each entry is the value of orthogonal projection of the pupil onto the basis set indexed by 'a*. At step S512, 
the weights of each respective discrete point in the illuminator profile from step S504 are combined with the corre- 
sponding orthogonal projection on the basis set as calculated in any one of steps SSIOjero, SSIO^pii, and S510(||tf> 
The solution to the combination perfonDed at step S512 renders the kemel A[ij] at step S514. 

15 [0083] The imaging of a simulation system In accordance with the exemplary embodiment of the present invention 
will now be described with respect to coherent imaging, incoherent imaging with the Hopkins model, and with respect 
to the basis set decomposition. 
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COHERENT IMAGING 
Pupil Plane 



[0084] The field at the mask exit plane <|>'"'*is given as the modulation of the incident field <|»'"-^ by the mask transfer 
matrix T, For an Inf inttesimally thin mask \a/il Tie effectively diagonal (For a basis that satisfiesFff i -i)> = (nlr) a diagonal 
25 matrix <rl Tlr) implies that <nl 77n*>=^il-'n*)): 

30 The field at the pupil entrance plane follows via equatk>n (39) as 

<<|>f;{ls^>=/ko<<|>[^flG^ls^> (47) 



= iko<*[s/lTG%Sp> (48) 



Substrate Plane 



40 [0085] The propagated image at the substrate plane follows according to equation (44) as the field propagated 
through vacuum segment 1 from the pupil-plane exit field 

1*1 



<*|^'lr^> - iko{^^^\GhrJ (49) 

[0086] The pupil-plane exit field itself is the pupil plane entrance field ^'[^ propagated by the lens transfer matrix 
</r 1 LM /(*). The substrate plane field thus follows as 

50 <4>ls)«0 = »Ko<4>r;j >t.^G^ Ir^) (50) 

[0087] The lens transfer matrix Is assunr^ed to be diagonal In Fourier space, (MU I A')= (k\U I A>8(*f,/f). Using equation 
(47) Ihe substrate plane field may be expressed as a function of the field at the mask input plane as 



^[4K = -k^<c|)S'lTG\^G^Ir^> (51) 
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Source Plane 

[0088] Coherent Imaging assumes a plane wave incident at the mask at Incidence angle d^. Illumination at this angle 
Is equivalent to Illumination at nomnai incidence with a phase modulation of the mask plane amplitude by /c^ ~ ^ sin 
'dg = koSg. The scalar field at the mask exit plane for arbitrary Incidence angle can thus be expressed In temris of normal 
Incidence as: 



10 



(52) 



IS 



20 



[0089] Under Kohler illumination the field (^^\ /) reflects a plane wave of nonnal incidence. The innaginary part of 
such a plane wave is a globally unifonn phase factor and may be factored out. What remains is the real part of <(>^;'l 
/) reflecting the weight of the originating source point in the illumination aperture: 



[01 



(C'l /)=iv(^3). 



?[0]* '/^'^yvsh (53) 
[0090] Using the above Illumination, the substrate plane field for point source located at ^^then follows as: 



<*n'^w>=-w(^s)»<o /dr^e ^ '"(r^lTG^L^G^r^) 



(54) 



25 



Fraunhofer approximation 

[0091] The above model expresses a Fresnel approximation of the lens system and is nonlinear. Several approxi- 
mations follow to linearize the above imaging system. Expanding equation (54) into its spectral components yields 



30 



35 



(4n) dod^ 



(55) 



(56) 
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X exp 



2d„ 



2f. 2d. [ 



(57) 



[0092] The argument of the exponential may be evaluated via: 



•45 



50 



In A = ikoSpi 



(58) 



[0093] In the last term, the sub-term associated with ^ will fall out when squaring the field to obtain the Intensity 
of the field at the substrate plane. Under the assumption ^fit the image at the substrate plane approximately mirrors 
the image at the mask plane, r^ - -Mr„, the sub-temn ^ can be eliminated for the same argument. 

DefpcusTerrn 



55 



[0094] The centerterm In the above described exponential relates to the focus condition of the system. For a perfectly 
focused system, this term disappears according to equation (34). For a system under small defocus of magnitude ^, 
a Taylor expansion yields: 
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1 1 



d, +? d, d^ 



(59) 



wherein a phase term accounting for defocus thus follows as: 



10 



(60) 



[0095] A more accurate defocus model follows by calculating the true optical path difference ("OPD'O of the Gaussian 
IS reference sphere: 



20 



(i>'^(c)=koca-<T.) 



(61) 
(62) 



Rest of Terms 

[0096] The first term in equation (58) contains the main imaging information and corresponds to a Fourier transfor- 
mation of the system. Introducing the scaled coordinates equation (35) this temi follows as: 



30 



35 



40 



is ( r -r ) 



(63) 

(64) 
(65) 



[0097] Defining 



45 



P(Sp):=P'(sp)exp 



so 



2d,' 



<; = P'(s,)exp 



2k, 



(66) 



the lens system equation (55) then follows as 



55 
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(67) 



Define <rjflr^:=<r„mr^4plAsp>:=<SplPls^. then 

(0'^«)=^('5'.)(s,|tP|r,) (68) 
INCOHERENT IMAGING: HOPKINS MODEL 

[0098] In the following, the "hat" notation is dropped and scaled coordinates according to equation (35) are assumed 
throughout. For incoherent imaging the Intensity at the substrate plane follows as a weighted summation over all In- 
dependent point sources ^ : 



I(r,)=Jdk,|((i>;;;|r«)|^ (69) 
= K(^.h''l^)'«'(^)'{^J7-p|r^) (70) 
= I ^l*'X*' I ^l^)'*''^x>''(^ \T\k){k\P\rJ)dkdk' (71) 



Let 



t(k'-k) = {k'\T\k) q' = fr'-«, <7 = *-s, (72) 



then 



l(r^, = Jdqdq-dk, <r„|P|q'+gs>t(q')w(S,)2l*(q)<q+g|P| r^> (73) 

= !dqdq't(q')t'(q)ldk/rJP\<f^.^^w(^^^(q4.sJiP\r„) (74) 
since <;ciPlft> is diagonal define P(ki :=<MP/ik>: 

I(r, ) = Jdq dq't(q')t *(q) Jdk, (r, I q' + )P * (q' + Ss Hh f (q + s J )P(q + s, ) (75) 

= idqdq^t(q')t'(q){r^ \q- .q)idk ^P'{q' ^s^]v{s^y p{q +s^) (76) 
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Define the Tratsmisston Cross Coeffifcfenf (TCC) function T{q',q) 

T(q'.q)= /dk,P*(q'+S,)w(S,)*P(q+S,) (77) 



5 

then 



10 jdqdq'{r^ | q') t ') C(<7'.^) «' W) | r^) (78) 

Assume a factorization exists sucli that 

15 

T(q'.q) = 2;;7;(qOA,nj(q) (79) 
where tiXq) forms an orthogonal basis. Combine equation (77) and equation (79): 

2;Tl*(q')Agni(q)=Jdk.P*(q' + §.)w(s.)^P(q + s.) (80) 



2S 
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[0099] Let Nf = (ti/ti;) be L^-nom of T|yw Left- and right multiply equation (80) with the orthogonal basis ti/y followed 
by an integration over the arguments: 

^=nJn/ "^"^ d«^'dMi(q')P*(q'->-S3)w (s^)^P(q+S3)Ti;(q) (81) 

[0100] An orthogonal pupil proiec^on coefficient ^Hk^ may then be defined as: 

Pi(K) :=jq-/dqP(q+C)n;(5) (82) 



l/dg'P(5*)n;(5* - h (83) 



[0101] In accordance with the present Invention, pXQ may be pre-caiculated for any possible pupil configuration, i. 
e.. considering non-limiting parameters such as NA, focus, aberration, etc. The radial part of the P(k) is typically a 
circular aperture. The phase part describes aberrations of the tens field. 

[0102] Using the orthogonal pupil projection, the coefficient matrix, or l<ernel, A»of equation (82) follows as: 



so 



A,j = /dkw(l<)p;(k)pj(lc) (84) 



Note that A/y is IHermitean. 



[0103] Stop S506 may be porfomried In accordance with above-described equations (20) and (1 9), and as illustrated 
with respect to FIG. 7. In a first step S702, the column a = 0 Is initialized in accordance with equations (11) or (20), 
55 wherein each entry In the column is seeded with an incomplete gamma function. In the next step S704, the values of 
column a =r 0 are propagated to the column to the right in accordance with equations (10) or (19), wherein each entry 
in the column is filled according to recurrence. After repeating step S704 until row m = 0 is full, the resulting indefinite 
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integrals, values for J[a](/f/) or J[al(x^, may be determined. As such, the present Invention generates a table of Indefinite 
Integrals based on a recun-ence over a seed array of Incomplete gamma functions. 
[0104] Steps S508 - S514 will now be described in more detail. 

^ Algorithms: Aerial imag e 

a. Pupil 

[01 05] Combining equation (82) and equation (66), the orthogonal pupil projection coefficient follows as 



10 
15 



P.(s.)=^/dqP'(q)exp 



ni(4-s.) (85) 



Approximation 1: Alserration IModel 



[0106J An aberration model is obtained by modulaUng the phase distortion over the pupil opening. Normalize the 
pupil coordinates with respect to NAs=1: 

A 

(86) 

25 

using Z4 = 2 ^- 1 the phase error K caused by defocus can be expressed as: 

• = |^«=7ko€'(Z4+2,) (87) 

• [01 07] The constant phase tenn 2^ may be ignored. Higher order temns of the phase error may be loaded extemalhf 
from a lens aberration map via a coefficient set a„. (a>„ := a^^. Define 

35 



'-Y<t>, (88) 
[0108] The orthogonal pupil projection coefficient then follows by numerical Integration of 

P.(ss,s) exp [i<J>(f,q)]ti;(q-sJ (89) 

[0109] For tenns that have only radial dependence (q :=lql) the integration may be simplified as: 

P3(Ss><;) =^/o"qdqexp[lO(f,;)]{'V^ii:(q-sJ (90) 

S5 Approximation 2: Zero Defocus 

[01 10] For lenses without aben-atlons, the pupil function is one over the domain of the aperture and zero outside 
The pupil aperture may then be expressed as a finite circular integration: 



so 
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[01 11 ] For a lens system that Is In-focus (C = 0) equation (85) reduces to a convolution of the pupit aperture with the 



spectral basis and p^ (q^.s g) 



10 



P.iO (92) 

t5 [0112] The circular aperture may be approximated by a superposition of rectangular domains S. as illustrated In FIG. 
3. FIQ. 10, on the other hand, Illustrates the improved method of the present invention, wherein for each domain B, 
the pupil convolution evaluates to 

pf(s>pf;(a,.,)p;;(S,,) os) 



25 
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[01 1 31 An efficient method for evaluating the above one-dimenstonal integral, will now be described with respect to 
FIG. 11. 

[01 1 41 From equation (1 7) the zero defocus case may be expressed as 

[01 1 51 Sampling over a finite grid a finite approximation to Afj follows as: 

A,. = A-XZ^ij/^;(^)p>,) (96) 

[01161 Note that (k), and thus p^Cs^), are purely real for even indices and purely Imaginary for odd Indices. The 
image amplitude Is evaluated using only the temis 2Re {A/j}. Thus, for In-focus conditions all entries In A/y that are 
purely Imaginary may be disregarded, as illustrated in FIG. 11 A, as sub-block 'D'. The remaining diagonal A and sub- 
blocks B and C may be organized in software as illustrated in FIG. 11B. 

Approximation 3: First Order Defocus 

[0117] A first order approximation towards defocus aberration follows by calculating a series expansion of Pa(s^ 
around ^ - 0 up to the first order tenn: 

[01181 The first term is given in equation (95). The second tenm may be evaluated as: 
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n;.(q-s.) (98) 



[01 19] For {; = 0 this term follows as: 



5«;p.(s..(; = 0) = ^ Jdqg'nKq-s.) (99) 
[0120] A first and second order correction to the image kernel follow as 



K = ^^s,i:-(3j[p:(s.)a,p,(sj+(a,p;(sj)p,(s.)] (100) 



A«> = A^s,i»v(sja^p;(iJa^p.(fJ (loi) 



[0121] Note that A(^)fj Is Hermitean with even-even and odd-odd elements purely imaginary entries while >4<2>^^ \q 
Hemiitean with even-odd elements purely imaginary. Calculation of the Image intensity yields a bilinear fomi with a 
coherent vector g. For binary, as well as ISO** phase shift masks, all entries G pare real. For /KD^thus only the 
even-odd and odd-even elements are of relevance while for A(^)g only the even-even and odd-odd elements contribute 
to the image intensity. 

[01 22J The first order con-ection to the Image Intensity then takes the form: 



l(r,0 . 2<'-|^>„.(^. lr)^(r|g„)A;';^>(g^lr>C (102) 



[0123] The first order defocus model can be utilized to rapidly calculate a number of Image planes throughout the 
thickness of the resist profile. As long as the focus deviation over the resist thickness is small (comparing the phase 
distortion In the outer pupil ring vs 2rt), the first order approximation gives a sufficiently accurate model of the defocus 
behavior through a resist film of finite thickness. 

Algorithm 4: Improved OPD integration scheme 

[0124] For a general solution of equation (91 ), the following steps may be used: 

1 . DIscretIze the domain of the pupil over a finite grid with indices k„ s„ 

2. Build a matrix of OPD values in the pupil plane M^^k^ky^Q := c/rc(/c)exp[/5>(/f^,/f^0] 

3. Build vector of basis function values In first dimension: A{Bi^, q^) := t| (q^) 

4. Calculate the following vector: 



9? = -——{(^qq^ej^p 

'^'^O^i o 



2k 



0 . 
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C (k,,ay)[sJ=^Mp(k,,ky,^)A(kj,-Sy,ay) 

k 

y 

(103) 

5. Use the above matrix of vectors to combine all calculations In the second dimension: 
Pa is ,s .fl )- — zAk -s ,c )c ik ,a % ] 

'if 

(104) 

incoherent Imaging: Hopiclns IHodei witli Diffusion 

[0125] In the discussion Immediately below, the Hopkins model is modified to Include the effects of photoactive 
^0 compound diffusion In the resist. An isotropic Gaussian diffusion model is assumed. Diffusion Is modeled as convolution 
with a Gaussian Icernel: 

25 <r\^\n :=-Ue (106) 

2no 

with 

30 

u = j2Di (106) 

[01 261 Assuming the photoactive compound concentration Is linear in /(/), diffusion may be approximated by replacing 
equation (73) with: 

35 

<''« '-^<'a) (107) 
[0127] Comparing the above to equation (70), this is equivalent to substitute 

40 

P'-^ZP (108) 
[0128] The intensity In equation (74) then follows as: 

45 

l(r^) = / dqdq't(q')t*(q) Jdk3(r^|ZP|q'+§3)w(23)2<q+gjPZ|r„> (109) 
[0129] The bracketed temn inside evaluates to: 

so 
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(rJZPl^' + i,) = ldkdk\r^\k){k\^k){k\P\q'^s^^^^ (110) 

5 

= \dk{rj^k){k\z\q'^s,)p{q'^s,) (111) 

/o [0130] If ( /y2l r) = /^r- 0. then 2 is diagonal In k-space. Thus: 

<rJZP|(7Va,> = <rJq'+S,>Z(q'+a^)P(g'+g^) (112) 

[01 31 ] Therefore, diffusion can be modeled by including the Fourier transf onu of the Gaussian kernel equation (1 05) 
into the calculation of the orthogonal pupil projection coefficient (c.f. equation (82)): 



= fJ-JdqZ(q-i-^)P( q4-K)n; ( q) (113) 



20 

[01 32] As relating to FIQ. 6: 
Point Decomposition 

25 

[0133] Define: 

(r^lgj) = /dk<r^ lk){tlk><klTij) (1 14) 

30 

[0134] Then equation (78) can be expressed as: 



35 



liO = ^M''^\9i){9j\r^) (115) 
(fif/l amounts to a convolution integral in real space: 

40 

{9i\ i) =fdr^}f\ r-0(flO* (116) 
[0135] Decompose the mask transmission function into a disjoint set 

45 

{t\r) = Z'pKi'') (117) 

P 

[01 36] Where tp is the complex transmission of the layout structure indicated by p and 1 ^ an Indicator function of unit 
value over the support of p. For Manhattan geometry's 1^ is rectangular and may be factored as 



V= V '1 p. (118) 



55 

and thus 
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{9<\r) = Jit;fdr'{7j,\r-r')l^{r) (119) 

P 



Reduction off Fieid Integral 



[01 37] The above expressions reduce the mask transmission function to a summation over an Indicator function 1 p. 
Furthermore, tor Manhattan geometry's, this Indicator function factors into the cartesian product 1p= Ip;^* 1 py. The 
10 next step is to simplify the Integration In equation (1 1 7) utilizing the above reduction of the mask transmission function. 
[01 38] Factorization of the indicator function is matched by the choice of the spectral representation as discussed 
above with respect to basis set decomposition. The basis t|/ Itself also factors into a Cartesian product ^„-r\^ 't]^ 
• This penults decomposition of the two-dimensional field integral of equation (117) Into a product of two one^mei?^ 
sional field integrals. Define: 



I„{dp,\dp,') := Jdx'(Ti:;jx-x')(l,Jx') (120) 

= ft; '^'R 1^-4 021) 



[0139] Where 

30 "^^x •= ^ ' Hx,o 



dp° := X - Q lower boundary (1 23) 



dpi := x-Pjj ^ upper boundary (124) 

35 [0140] Using the Gamma-function factorization as discussed above with respect to the spectral basis, the above 
integration may be reduced to a recurrence over the seml-lndefinite integrals JJi)^ via equations (4), (9) and (10). 
[0141] The above expression equation (11 7) then simplifies to: 



(^-i'-) = E';«..(*'x.4pi)-i„,(4p;.«5p;) (125) 



[0142] A particular advantage of this representation follows when several layout features line up along a Cartesian 
^5 dimension. Assume for example q e to be a subset of features that are vertically stacked up and have identical 
transmission t^y. Then for alt gthe values of q^Q = and q^^^ ^ are constant. The above expression then reduces to: 



(Pf'lr) = t;,,l„^(dc^dc■)X^n,(clq^dq'J (126) 



Intensity Calibration 

55 [0143] For afully transparent mask, I,p^p{i)= 1, =1 and equation (117) yields: 
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{qT^ d==idr'{t\, \ r-r*) (127) 

= ^,(-*^.«)ln3,(-~.~) (128) 



Dark Field Reticle (binary) 



[0144] A binary darkfleld mask has no background transmission while features have a transmission of 1 . For Man- 
^0 hattan geometry features factor according to: 



1p=1p^»V^ (129) 



C lear Field Reticle (binary) 
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[01 451 For a binary clear field mask features have zero transmission (OJ while the background has unit transmission. 
Note that 



E^P-^S^P' ^ 1 (130) 



[0146] Thus: 



igf^l^ = {g7 l^<9l^l^ (131) 

30 

Clear Field Mask (binary attenuated phase shift) 

[01 47] For an attenuated phase shift mask with 1 80° phase sh rft and attenuation factor a, the transmission of features 
is ( a) while the l)ackground transmission is +1 . Thus 



35 



40 



('I'-) = Zi;('-)-'^Zi;- (132) 



[01 48] Note again that 1 +p and 1 completely cover the simulation domain: 



= 1 (133) 

45 p p' 

[0149] Thus 

(^i'') = (134) 

= l-(^ + l)Z^P' (135) 



so 



55 
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Clear Field Reticle (ternary attenuated phase shift) 

[01 50] A ternary attPSM clear field reticle consists of three transmission domains: a 0** phase fully transmisslve clear 
field region 1^, a 180** phase shifted region 1 ' with transmission a. and a chrome blocking region 0- A full partition of 
5 the mask area follows by the disjoint union of all three domains: 



= 1 (136) 



10 ^ ^ , 

[01 51 ] The mask transmission function can be expressed as 



('I'-) = HKirhaJ^r, (137) 
p tf 

2o = J-Zl^-'^Sl^-So, (138) 

= \-ia^X)YXp-Y.^, (139) 

25 P' 9 

[01 52] The coherent field vector follows as 

30 {qT"^^ Ia) = {g; l/)-(a 4- 1)<pf^ Ir- {gf^^ li) (140) 

Dark Field Reticle (ternary attenuated phase shift) 

[01 53] The domain partition of a ternary attPSM darkfield mask is identical to the domain partition of a ternary attPSM 
55 Clearfield mask. The mask transmission function however is expressed In terms of 1+^ and 1"^: 

('I'-) = (141) 

40 P p' 

and thus the coherent field vectors follows as: 

45 ^gdanpsr^c ,^ ^ ,^ ^^^^^ 

[01 54J For an attPS!\/l contact mask, the number of rectangles can be reduced through a different linear combination 
of 1 + and 1 For a contact at index p define an 'outer frame region' 1/ that includes both the 0** and the 1 SO" region as: 

so 

= (143) 

the mask transmission function then follows as 

55 
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(144) 



P 



(145) 



P 



(146) 



P 



the coherent field vector can then be expressed as 



(flff-'*^l/) = (Ua)<p;i/).a<pfl,> 



(147) 



Fast Implementation of Zero and OPD Model 

[01551 The present Invention provides an efficient method for calculating orthogonal pupil projections in the zero- 
defocused models, as Illustrated In FIQ. 10. Equation (92) is used to calculate an integral over the circular pupil aperture 
in (x,y) space. As Illustrated in FIG. 10, the center of the aperture is shifted by S ^. The method for performing the 
algorithm as illustrated in FIG. 10 is as foBows. First is an angular step-size of 45^ wherein the point at the circle arc 
intersecting a ray emanating from the center Is located. Next is to find the three mirror points, across the center and 
across the x and y axis through the center of the circle. These points are connected to form the rectangle Bq. Then the 
contribution of to p[a] according to the equations (95) and (17) is determined. 

[0156J The following steps are then performed until the resulting rectangle becomes smaller than a discrete Increment 
of the pupil grid. The angular step-size is reduced by half. Again, the point at the circle arc Intersecting a ray emanating 
from the center at the cun^ent angle is detenmined. Again the symmetrical point is detennined and the new set of 
rectangles is detemnined. Then the contribution of the new rectangles to pfa] according to the equations (95) and (17) 

is detennined. 

[0157] The method of calculating the area within the pupil as described with reference to FIG. 10. is much more 
efficient than the prior art system of FIG. 3. This is because In the method of FIG. 10. every discrete portion of the pupil 
is not calculated, thus a large amount of calculation time is saved. Specifically, the method of FIG. 10 uses a smaller 
number of rectangles to approximate a circular pupil shape. 

[01 581 Step S408 will now be further described with respect to FIG. 6. The Initialization step S602 comprises three 
steps, S604, S606. and S608. At step S604 a discrete grid of the mask is specified for the geometric sampling. Methods 
for perfomiing this step are known In the art. Step S606, however. Is perfonned contrary to that of the prior art. 
[0159] As opposed to an-anging the data of the geometry for discrete sampling points on the grid of the mask as In 
the prior art systems, the present invention tabulates a one-dimensional an-ay In other words, the prior art method 
uses a one-dimensional array for the spectral indices and a two-dimensional array for the geometrical Indices, thereby 
resulting in an array of two-dimensional tables. On the other hand, the method of the present invention uses a one- 
dimensional array for the spectral indices and a one-dimensional an-ay for the geometrical Indices, thereby resulting 
in an array of one-dimensional tables. 

[01601 Specifically, an array of seed-elements for the spectral basis of each discrete point within a mask grid corre- 
sponding to the sample area Is tabulated. Each seed-element Is an Incomplete gamma function J[a], such as Is de- 
scribed above with respect to equation (11) and equation (10). The final Initialization step S608 calculates the clear- 
field exposure reference and saves the value for subsequent exposure dose calibration with equation (128) and equa- 



[01 61 1 Step 861 0 will now be described in more detail. In order for an aerial image of a mask pattern to be simulated, 
the geometry of the respective mask must be decomposed into non-overiapping rectangles. As further illustrated In 
FIG. 8, item 800 is a prcdetennlned proximity window containing an Imago of Item 802, a sample mask comprising a 
plurality of sample features 804-808. Each feature is decomposed into non-overiappIng rectangles. For example, fea- 
ture 804 is decomposed Into non-overiapping rectangles 81 0 and 81 2. 

[0162] A circular proximity window of radius may be defined as the domain for which the normalized edge 
variation of the zero-order basis <gol i) remains larger than or equal to G^^.^. The boundary of this domain follows as: 
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<tlSlr=0) 



Thus 



1, prox 



(149) 



[0163] In essence: the proximity window defines the set of alt rectangles that affect the sannpled area. Specifically, 
In this example, as illustrated in FIG. 8, features 814-818 do not have an affect on the aerial image of features 804-808. 
Therefore, features 814-818 are not within the proximity window. 

[0164] Step 861 2 will now be described in more detail. Once the features within the sampling window have been 
IS decomposed into non-overlapping rectangles, the corresponding orthogonal projections of the features of the mask 
onto the basis set are calculated in accordance with equation (125) and equatbn (9). 

[01 65] Step S61 4 will now be descrit>ed in more detail. Once the orthogonal projections have been calculated In step 
S612, corrections may need to be applied, for example, depending on the mask. If the mask is a clear field mask, then 
equation (131) may be used lo correct the orthogonal projections. Further, If the mask is an attenuated phase-shift 
20 mask, then equation (1 40) may be used to correct the orthogonal projections. Finally, if the mask Is a ternary attenuated 
phase-shift mask, then equation (142) may be used to correct the orthogonal projections. Once the correctk>ns have 
been applied In step S614, the final corrected orthogonal projection of the mask transmission function g[4 is attained. 



Intensity Sunrwnatlon 



[0166] Returning back to FIG. 4, once steps S404 and S408 have been performed, the image intensity distribution 
may be calculated at step S410. Step S410 will now be further explained with respect to FIG. 9. As illustrated In FIG. 
9, at step S906, the A[ij] from step S902 is multiplied with. the g[4 from step S904, in a bilinear form. This cateuiatton 
provides the image intensity /(#) at point r. Furthermore, the intensity gradient at ponit rmay be calculated at step S908. 
30 As Illustrated in step S912, the intensity gradient may be measured with respect to the gradient In either the x or y 
coordinate axis. 

[01 67] Specif toally. the Image intensity /(/) at a particular sampling point, such as In step 8906, can then be calculated 



= Z{'-'9M^{9„\r) (150) 



a./(r) = 2ReX^^(r|g„)(:^(gj„,.,Jr)-(g,„,„,|r)) (151) 

nm 

45 [0168] A numerical Implementation may take advantage of the hennlticlty of A^. The basis r\/if) Is real. For binary 
masks, as well as 180° phase shift masks, tp is addilionaily real. Thus, the number of summation terms may be cut 



approximately In half: 



Dose Calibration 



[0169] For a computational implementation, a calibration may be used against a Clearfield area for absolute exposure 
calculations. For example, using equation (128), the intensity of the clear-field exposed area then follows according to 
equation (150) as: 
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(153) 



5 



[0170] Therefore, given an exposure dose D, the local intensity sampled at point roan then be calculated as: 



10 



■^W = -^i;^«<''ki><gj|r> 



(154) 



Equation (150) calculates A which is the Intensity calculated by the non-nonnalized algorithm for 100% exposure. 
IS Therefore, In order to calibrate the systems to give an exposure dose value D as specified by the user, a prefactor of 
D/ Is applied in step S906. 

[0171] The system and method of the present Invention may be Incorporated into many applications. Non-llmlting 
examples Include: core simulation engine as part of an Model OPC calculation; simulation engine for predicting and 
optimizing fabrication parameters, such as NA, illumination, etc., for drawn layout features; simulation engine for pre- 
20 dieting and optimizing geometric assist features, such as serlfe, scattering bars, etc, to enhance the printability and 
yield of drawn layout features; simulation engine to predict the printing of defects on a mask for a specified fabrication 
process; and simulation engine to predict and optimize the balance of lens aberration components to optimize the 
printability and yield of drawn layout features. 

[0172] Although specific reference may be made in the text to the use of lithographic projection apparatus in the 
2s manufacture of integrated circuits, it should be explicitly understood that such apparatus has many other possible 
applications. For example, it may be employed In the manufacture of integrated optical systems, guidanceand detection 
pattems for magnetic domain memories, liquid-crystal display panels, thin-film magnetic heads, etc. The skilled artisan 
will appreciated that, in the context of such alternative applications, any use of the temns "reticle" or "wafer" In this text 
should be considered as being replaced by the more general terms "mask" or "substrate", respectively. 
30 [0173] Furthermore, although this text has concentrated on lithographic apparatus and methods whereby a mask is 
used to pattern the radiation beam entering the projection system, It should be noted that the invention presented here 
should be seen In the broader context of lithographic apparatus and methods employing generic "patterning means" 
to pattern the said radiation beam. The temn "patterning means" as here employed refers broadly to means that can 
be used to endow an incoming radiation beam with a patterned cross-section, corresponding to a pattem that Is to be 
3S created in a target portion of the substrate; the temi "light vah^e" has also been used In this context. Generally, the said 
pattem will correspond to a particular functional layer in a devfce being created in the target portion, such as an 39 
integrated circuit or other device. Besides a mask (whether transmlssive or reflective) on a mask table, such patterning 
means Include the following exemplary embodiments: 

^0 (1 ) A programmable mln-or an^ay. One exannple of such a device Is a matrix-addressable surface having a vlscoe- 

lastic control layer and a reflective surface. The basic principle behind such an apparatus is that (for example) 
addressed areas of the reflective surface reflect Incident light as diffracted light, whereas unaddressed areas reflect 
Incident light as undiffracted light. Using an appropriate filter, the said undiffracted light can be filtered out of the 
reflected beam, leaving only the diffracted light behind: in this manner, the beam becomes patterned according to 
the addressing pattem of the matrlx-adressable surface. An alternative embodiment of a programmable mirror 
array employs a matrix an-angement of tiny mirrors, each of which can be Individually tflted about an axis by applying 
a suitable localized electric field or by employing piezoelectric means. Once again, the mirrors are matrix-address- 
able, such that addressed mirrors will reflect an incoming radiation beam In a different direction to unaddressed 
mirrors; in this manner, the reflected beam is pattemed according to the addressing pattern of the matrix-adressable 

so mirrors. The required matrix addressing can be perfomned using suitable electronic/computer means. More Infor- 

mation on mirror arrays as here referred to can be gleaned, for example, from United States Patents US 5,296,891 
and US 5,523,1 93, and PCT patent applications WO 98/38597 and WO 98/33096, which are Incorporated herein 
by reference. 

55 (2) A programmable LCD array. An example of such a construction Is given in United States Patent US 6,229,872, 

which Is incorporated herein by reference. 

[0174] Any reference to the term "mask" in the Claims and/or Description should be Interpreted as encompassing 



27 



EP 1202 119 A1 

the term "patterning means" as hereabove descHbed. 

[0175J In summary, although certain specific embodiments of the present invention have been disclosed, the present 
Invention may be embodied in other forms without departing from the spirit or essential characteristics thereof. The 
present embodiments are therefor to be considered In ail respects as illustrative and not restrictive, the scope of the 
invention being indicated by the appended claims, and alt changes that come within the meaning and range of equiv- 
alency of the claims are therefore intended to be embraced therein. 



Claims 



10 



1. A method for simulating an aerial image projected from an optical system, said optical system including a pupil 
and a masl( plane, the method comprising the steps of: 

providing a mask to the mask plane; 
obtaining parameters for said optical system; 

calculating a kernel based on an orthogonal pupil projection of said parameters of said optical system onto a 
basis set; 

obtaining parameters of said mask; 

calculating a vector based on an orthogonal mask projeclion of said parameters of said mask onto a basis set; 
^ calculating a field Intensity distribution using said kernel and said vector, and 

obtaining aerial image data from said field intensity distribution. 

2. The method of claim 1 , wherein said parameters for said optical system Include aberrations. 

3. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
optica] system includes the step of tabulating an array of incomplete gamma functions corresponding to respective 
points in said pupil of said optical system. 

4. The method of claim 1 or 2, wherein saki step of cak^ulating a kernel corresponding to said parameters of said 
30 optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in said pupil of said optical system, 

wherein said optical system is In-focus. 

5. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
35 optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients con«- 

sponding to respective points in said pupil of said optical system, 

wherein said optical system Is either not In-focus or has aberrations. 

6. The method of claim 1 or 2, wherein said step of calculating a kemel corresponding to said parameters of said 
-^0 optical system further includes the step of tabulating an array of orthogonal pupil projection coeffteients corre- 
sponding to respective points In said pupil of said optical system, 

wherein said optical system accounts for effects of photoactive compound diffusion in a resist layer present 
on a substrate onto which an image of the mask Is to be projected. 

-'^ 7. The method of claim 4, 5 or 6, further comprising the step of combining sample weights of an illuminator profile of 
said optical system with said array of orthogonal pupil projection coefficients. 

8. The method of claim 1 , wherein said step of calculating a vector corresponding to said parameters of said msisk 
further includes the step of specifying a proximity window within said mask for geometric sampling. 

50 

9. The method of claim 1 , wherein said step of calculating a vector corresponding to said parameters of said mask 
further includes the step of decomposing a geometric pattern of said mask Into a disjoint set of rectangles and 
tabulating an array of projections of the rectangles within a proximity window. 

55 10. The method of claim 9, wherein said step of calculating a vector con^esponding to said parameters of said mask 
further includes the step of con-ecting the array of projections of the rectangles based on the type of mask. 

11. A method for simulating an aerial image projected from an optical system, said optical system including a pupil 
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and a mask plane, the method comprising the steps of: 

providing a mask to the mask plane; 
obtaining parameters for said optical system; 
obtaining parameters of said mask; 

orthogonally projecting polynomials, said polynomials Including polynomials associated with said parameters 
of said optical system, and 

approximating the transmission cross-congelation function associated with said optical system based on said 
orthogonal projection of polynomials. 

12. A processing device operable to simulate an aerial Image projected from an optical system, said optical system 
including a pupil and a mask plane, said processing device comprising: 

a first parameter obtainer for obtaining parameters of said optical system; 

a first calculator for calculating a kernel based on an orthogonal pupil projection of said parameters of said 
optteal system onto a basis set; 

a second parameter obtainer for obtaining parameters of a mask provided to said mask plane; 

a second calculator for calculating a vector based on an orthogonal mask projection of said parameters of 

said mask onto a basis set; 

a third calculator for cak^ulating a field intensity distributton using said kernel and said vector, and 
an aerial image obtainer for obtaining aerial image data from said field intensity distribution. 

13. The processing device of claim 12, wherein the first, second, and third calculatore are the same cateulator. 

-14. A processing device operable to simulate an aerial image projected from an optical system, said optical system 
including a pupil and a mask plane, said processing device comprising: 

a mask within said mask plane; 

a first parameter obtainer for obtaining parameters of said optical system; 
a second parameter obtainer for obtaining parameters of said mask; 

a first calculator for orthogonally projecting polynomials, said polynomials including polynomials associated 
with said parameters of said optical system and polynomials associated with parameters of said optical system 
and 

a second calculator for approximating the transmission cross-correlation function associated with said optical 
system based on said orthogonal projection of polynomials. 

15. A device manufacturing method comprising the steps of: 

(a) providing a substrate that is at least partially covered by a layer of radiation-sensitive material: 

(b) providing a mask that contains a pattern; 

(c) using a projection beam of radiation and an optical system to project an image of at least part of the mask 
onto a target portion of the layer of radiation-sensitive material, 

whereby, prior to performing step (c), an aerial image to be projected from said optical system is simulated 
using a method comprising the steps of: 

obtaining parameters for said optical system, which comprises a pupil and a mask plane; 

calculating a kernel based on an orthogonal pupil projection of said parameters of said optical system 

onto a basis set; 

obtaining parameters of said mask, provided at said mask plane; 

cateulating a vector based on an orthogonal mask projection of said parameters of said mask onto a basis 
set; 

calculating a field intensity distribution using said kernel and said vector, and 
obtaining aerial image data from said field intensity distribution. 

16. A method of creating a data set representing a mask image, the method comprising the steps of: 

providing a data set representing a proposed mask image; 

simulating the aerial image of said proposed mask Image using the method of any one of claims 1 to 11 or the 
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processing device of any one of clainns 12 to 14; and 

adjusting or verifying said data set using the results of said step of simulating. 

17. A method of manufacturing a masic comprising manufacturing a mask using a data set created by the method of 
claim 16. 

1 8. A computer program comprising code means for causing a computer to operate to perfomi the steps of the method 
of any one of claims 1 to 1 1 . 
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Fig.2. 
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Fig.3. 
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Fig.9. 
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